On expansion of equal-time relativistic two-body wave equations in powers of 1/c 

to higher orders 
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Based on an extension of the Foldy-Wouthuysen method to two-body equations, the problem 
of expansion of equal-time relativistic equations for two Dirac particles in powers of 1/c to higher 
orders is considered. For the case of two particles with unequal masses, the transformed Hamiltonian 
in a general even-even form is obtained to order 1/c 4 . It is found that certain extra terms, which 
can be eliminated by an additional unitary transformation, arise in the expansion in higher orders, 
depending on the order of application of the generating functions in the first iteration. As examples 
for illustration, the Breit equation and the Salpeter equation with the Breit interaction are taken 
and their reduction to approximate forms including all the l/c 4 -order terms is carried out using the 
method under consideration. The obtained results may be applied for the nonrelativistic expansion 
of two-body wave equations with various interaction potentials to higher orders, for the investigation 
of their features and symmetries, and may also be useful in the study of light atoms. 

PACS numbers: 02.30.Mv, 03.65.Pm, 31.15.Md, 31.30.Jv 



I. INTRODUCTION 

In spite of considerable achievements in quantum field 
theory, particularly in quantum electrodynamics (QED), 
relativistic and semirelativistic two-body wave equations 
are widely used to learn about various relativistic effects 
in quantum systems and also to study and calculate their 
energy spectra (see, e.g., pj). In many cases, however, to 
obtain an acceptable description of a two-body system, it 
is not essential to solve the original equations of motion 
and one can restrict oneself to the consideration of their 
approximate forms with a small parameter, such as the 
inverse speed of light. As a rule, for many applications it 
is sufficient to restrict oneself to an expansion of the wave 
equations in powers of 1/c to second order. Thus, e.g., 
the well-known Breit-Fermi Hamiltonian and its special 
case for a vector exchange, which can be derived from the 
nonrelativistic reduction of different two-body equations, 
are often used to describe spectra of atomic, nuclear, and 
quark systems [HI]. 

A relatively simple description of a system of two spin- 
half particles, including relativistic effects and its energy 
spectrum, can be given with the help of equations of the 
Breit type. In order to expand the equations, the Foldy- 
Wouthuysen (FW) method Q generalized by Chraplyvy 
to the two-body problem is usually applied [g-[lj. 

According to this method, one represents a relativistic 
two-body Hamiltonian in such a general form: 



elements, and can be written as direct products of four- 
by-four matrices of each particle. In order to remove the 
undesirable terms (even-odd, etc.), the Hamiltonian must 
be subjected to canonical transformations of the type 



e lS He~ lS 



= H + t[S,H] 



2! 



[S, [S,H]} 



^-[S, [S, [S,H]}} 



(2) 



provided the generating functions S are suitably-chosen 
Hermitian operators, and where the inverse speed of light 
(or, equivalently, the inverse masses) can be taken as an 
expansion parameter. 

Though in the case of interacting particles one cannot 
determine a unitary transformation that cancels all the 
undesirable terms in all orders in principle, the original 
Hamiltonian can be converted by the procedure ^ into 
an even-even operator to any desired accuracy. There are 
a lot of sets of S @, @| to reduce a two-body Hamiltonian 
to approximate forms. In this paper, we use the genera- 
tors in a simple form, which, for the first iteration, read 
as follows: 
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-(SO), 



(3a) 



(3b) 



H = p l m 1 c 2 +hm 2 c 2 + (££) + (EO) + (0£) + (00). (1) 

This is a sum of the two "large" terms /3imic 2 + fizmic 1 
and even-even, even-odd, odd-even, and odd-odd terms, 
respectively. They are, in fact, matrices of 16 x 16 = 256 



'Electronic address: turovsky@bitp.kiev.ua 



Son — 



2(m\ — mfyc 2 



(CO). 



(3c) 



The generating functions for the next iterations differ 
from those in Eq. ((3]) by the structure of their odd-even, 
even-odd, and odd-odd factors. These generators enable 
one to reduce a sixteen-component two-body equation 
to its four-component approximate forms relative to the 
chosen energy states of the two particles. They were first 
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introduced by Chraplyvy in [6] to convert the two-body 
Hamiltonian ([1]) into an even-even operator to order 1 /c 2 
(under the assumption that (££), (00) are of order c°, 
and (0£), (£0) of order c 1 ), 



PlTTliC 2 + /?2W2C 2 + (££) 
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(4d) 
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(0£)(£0) 2 (0£) (4i) 



8mim2(mj — mfyc 4 
x [[(£0),(0£)],(00)] 



(4j) 



Because of the form of (|3c|) this prescription is applicable 
only for the case of unequal masses, however, it enables 
one to get four forms of the reduced Hamiltonian relative 
to the four different energy states of the system. 

Equation (|4]) was applied by Barker and Glover to the 
Coulomb case in order to transform the Hermitian part 
of the three-dimensional Bethe-Salpeter equation (the 
Salpeter equation [12[) and the Breit equation to sixteen- 
component approximate forms to order 1 /c 2 in [8| , where 
a comparative analysis of the transformed Hamiltonians 
for these equations was carried out and significant differ- 
ences between them were noted. 

Yet, a certain further expansion of the relativistic two- 
body wave equations to higher orders in 1/c may be of 
some interest as well, and not only from the theoretical 
or mathematical point of view. The improvement of ex- 



perimental techniques for the study of energy spectra in 
atomic systems, first of all, in hydrogen-like atoms, muo- 
nium, and positronium, allows one to make precise mea- 
surements of their energy levels [l3l - [l5| . They usually 
apply QED to learn about the a 5 mc 2 , a 6 mc 2 corrections 
to the energy (see, e.g., [l6|-[24j and references therein), 
however, the problem of derivation of the higher-order 
Hamiltonian of an arbitrary light atom remains difficult. 
Nevertheless, an expansion of relativistic two-body wave 
equations to order 1/c 4 provides a rather straightforward 
derivation of an effective a 6 mc 2 Hamiltonian, though it 
gives incomplete treatment of relativistic and quantum 
field effects. Still, it enables us to get some information 
about the structure of terms contributing to a 6 mc 2 for 
the energy eigenvalues. The nonrelativistic expansion to 
higher orders can also help to investigate the equations 
of motion themselves, their features and symmetries, to 
answer questions on relativistic effects they involve, what 
energies they work to, and so on. 

This paper focuses on the problem of expansion of the 
equal-time relativistic wave equations for two Dirac par- 
ticles to order 1 /c 4 by means of the extension of the FW 
method to two-body systems, using the generating func- 
tions ([3j) , and is organized as follows. Section [TTI deals 
with the two-body Hamiltonian transformed to higher 
orders, for which all the l/c 4 -order terms are found, and 
which is a continuation of the expansion Q in the case of 
commutation of (0£) and (£0). It occurs that the form 
of the higher-order part of the transformed Hamiltonian 
depends on the order of application of the functions (J3|), 
namely, it can involve certain extra terms having mass 
differences in the denominators even if both particles are 
in a positive energy state. Here, we are concerned with 
an additional unitary transformation canceling terms of 
this type as well. 

In Section IIII1 we work out all the terms of order 1/c 4 
in the expansion of the Breit equation and the Salpeter 
equation with the Breit interaction. In addition, we also 
discuss the possibility of modification of the l/c 4 -order 
part of the transformed Hamiltonians with the help of 
unitary transformations. Finally, Section IIVI contains a 
summary and main conclusions of this article. 



II. HAMILTONIAN TRANSFORMED TO 
HIGHER ORDERS 

Proceeding with the procedure of transformation of the 
Hamiltonian ([1}, using the generators ([3]) and with due 
regard for the commutation relation 



[(O£),(£O)} = 0, 



(•5) 



which causes a considerable simplification in £T tr , we get 
new even-even terms coming after the terms written in 
Eq. ([3]), and which are of lower order of magnitude. They 
form the higher-order, with respect to 1/c, part of the 
transformed Hamiltonian. 
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Thus one obtains such a prescription for the transformation of H: 



[( " } ' ( °' )]] + 8^ [( ' 0) ' mi {£0)]] + ^k? [{OE) > [{£ ° l {oo)] 
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■ iAw"-^ ' (0g> - KO£) ' [(0£) - (00)1+1 +1+ - S!?w-^ |w [(00) - <0£)21+l+ 



ftmi-ftm2 iv/n/n-> rcc/n-i r/r/m //oa\i i i ftlftrnj - ftm 2 ) 
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16TOi K - ™I) C6 K°°)' K £ °>' K fi °)' - i^- ^ K 00 ^ (f 0) ]+]+ (6d) 

MSO), (ff)], [(OS), (00)]+] + o Jl J [(OS), (££)], [(SO), (00)] + ] (6e) 



%m 1 m\& ' ' ' 8mfm 2 c 6 



384to 4 c 



384to*c 



^{[(Of ), [(OS), [(Of), [(Of), (ff )]]]] + 32[(0f ) 3 , [(Of), (ff )]]} 

Kji{K£0), [(SO), [(SO), [(SO), (ff )]]]] + 32[(f Of, [(SO), (ff )]]} (6g) 



' '(OS),[(OS),[(SO),[(SO),(SS)]]]\ (6h) 



647?i 2 m 2 c 8 

/3i/3 2 



)c8 -{[(^ ), [(Of), \(OS), [(SO), (00)]+]+]+]+ + 8[(Of ) 3 , [(fO), (00)]+]+} 
" " r {[(Of),[(fO),[(fO),[(fO),(00)]+]+]+]+ + 8[(Of),[(fO) 3 ,(00)]+]+} (6i) 



3„8 



96m 

" F(Of) 6 + T^Tn(£0) 6 . (6j) 



The FW method allows one to expand the Hamiltonian 
to any desired degree of approximation, keeping its Her- 
mitian character. Here, we assume that (ff ), (OO) are 
of order c°, and (Of), (f O) of order c 1 , and we retain 
the terms up to order 1/c 4 . So, under this assumption, 
expression (|6]) represents the transformed Hamiltonian 
approximate out to fourth order. It can be divided into 
two parts. The first part consists of the terms (|6a|) and 
is the transformed Hamiltonian to order 1/c 2 , into which 



the expression Q goes over under the commutation re- 
lation ([5]). The second part, which we will refer to as the 
higher-order transformed Hamiltonian, involves all the 
l/c 4 -order terms (|6b[..i). It is a sixteen-component equa- 
tion, and as usual, to obtain its four-component forms, 
i.e. reduced Hamiltonians, one has to put /3i = fa = ±1 
or f3\ = — fa = ±1. Note that if the inverse mass is con- 
sidered as an expansion parameter, the terms (|6bl ..e). 
which are nonlinear in (ff ), (OO), are of order 1/m . 
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Many of them have a mass difference in the denomina- 
tors. The terms ([BgJ h, i), which are linear in {££), (00) 7 
are of order 1/m 4 . The term (|6fj). which consists of the 
even-even and odd-odd operators, is the only one of order 
1 /to 2 in the higher-order transformed Hamiltonian. 

For the case when the mass of one particle becomes 
considerably great in compare with the mass of the other 
one, i.e. m-i — > oo or m 2 — ► 00, equation ([B]) goes over 
into the corresponding formula for a single Dirac particle 
in external fields: 



H tT = f3mc 2 +£ + J^0 2 



1 



P 



+ - 



8to 3 c 6 
1 



o- 



2mc 2 

13 



8m 2 c 4 
[0,£] 2 



[O, [£,0]} 



1 



8m 3 c 6 
O, [O, [O, [0,£]}}} 
P 



[0 s , [o,£]\ 



12to 4 c 8L ' l ' JJ 16m 5 c 10 
which is the transformation of the Hamiltonian 

H = I3mc 2 +£ + 0, 

with the use of the generator 



(7) 



2m& 



O. 



Under this condition only the terms of (|6b| g, j) remain in 
the higher-order transformed Hamiltonian and pass into 
the last four terms in Eq. ([7]); the rest of the commutators 
and anticommutators which involve the {OO) terms in 
the expression ^ vanish. 

The transformed Hamiltonian written out in Eq. Q 
to second order is derived with the use of the generating 
functions ^ and it is not very important which of them is 
used first in the first iteration in the series J2|); regardless 
of their order of application in this iteration, the same 
expression in the form of Q will be obtained. In other 
words, it is not important which of the undesirable terms 
in Eq. ([T]) will be destroyed first. But this statement is 
correct only if one needs to get the expansion up to the 
terms written in Eq. (|4]), or the same, to second order in 
1/c under our assumption. The use of S oe or S eo first in 
the sequence in the first iteration leads to the transformed 
Hamiltonian in the form ([5]), however, if one applies S 00 
first (instead of S oe or S eo ) in ([2]), certain "extra" terms 
of fourth order will arise in the transformed Hamiltonian 
in addition to those in Eq. ([BJ; namely, 



[[(0£)(£0),(00)},(££)} 



ftm 2 - fawi 
8raiTO2(m 2 — m^c 6 

8mim 2 (m 2 — m^c 6 



x[(£0){00){0£) - (0£)(00)(£0),{££)] 



m% - P\Pim\ 
16m 2 TO 2 (m 2 — m 2 )c s 

/3i/3 2 to 2 - mi 
16miTO 2 (TO 2 — m^c 8 

to 2 + ^ifemi 
16m 2 TO2(m 2 — to 2 ,)*: 8 



[[(0£)(£0),(00)],(0£) 2 ] 
[[(0£)(£0),(00)},(£0) 2 } 



x[(£0)(00)(0£) - (0£)(00)(£0), [0£f 



/3i/3 2 to 2 + mi 



16miTO 2 (TO 2 — m 2 )c 8 

x[(£0)(00)(0£) - (0£)(00)(£0), [£Of]. (8) 

As follows from this equation, the terms with mass dif- 
ferences in the denominators appear even though both 
particles are in positive or negative energy states (those 
correspond to setting j3\ — /3 2 = ±1)- Such terms, but 
with other numerical factors, also appear in H tr if one 
takes the sum S = S oe + S eo + S 00 as a generating func- 
tion in the series © . Obviously, some extra terms enter 
into the expansion and in sixth order and higher. This 
dependence of the form for the higher-order part of the 
transformed Hamiltonian on the order of application of 
the generators may be an important feature of the expan- 
sion to higher orders, and makes the difference between 
the second-order and higher-order expansions. 

Still, the extra terms can be eliminated given that 
the transformed Hamiltonian is subjected to an addi- 
tional unitary transformation like ^ with a generating 
function in the form of a Hermitian even-even operator, 
which we represent in such a manner: 



z(/3ito 2 - /3 2 TOi) 
8toito 2 (to 2 — to 2 )c 6 

ijjhmz + ftjTOi) 
8toito 2 (to 2 -to 2 )c 6 

x{(0£)(00)(£0) 



[(00),(0£)(£0)} (9a) 



{£0){00)(0£)j. 



(9b) 



Taking into account its higher order, one can retain only 
the first two terms in the series © , while the rest of the 
terms can be discarded because they are of lower order 
of magnitude, 



e ±/ tr e 



H tx +i[S ee ,H tI }. 



(10) 



Since Sqq I1&S an even-even form, it commutes with the 
two large terms from i? tr : 



[S ee , Pimic 2 + p 2 m 2 c 2 



0. 



(11) 



Actually, to remove the terms ©, it is quite convenient 
to retain the members of order c° in H tr standing in the 
commutator in (I10|) and to omit the rest. It is sufficient 
to get the terms which coincide with the ones in Eq. © 
up to a sign. One sees that the operators (|5a|) and (|9b[) 
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act separately. The former destroys the first, third, and 
fourth terms in Eq. ([5]), and the latter destroys the rest. 
We note that in general while removing the extra terms, 
this procedure gives rise to new ones instead, but all of 
them are of lower order of magnitude. 

Thus the generating function S ee allows us to modify 
the higher-order transformed Hamiltonian, subtracting 
(or removing) the terms (|5|). One can express it in terms 
of the original operators S oe , S eo , and S 00 in a convenient 
brief form: 



See — ['S'oei [<^e03 *~*oo]] 



(12) 



In conclusion of this section we should point out that 
if the commutation relation ([5]) had not been taken into 
account, the expression for H tr would be much lengthier 
than that in Eq. (J6j) and many other terms would appear 
in it as well. Moreover, the form of S ee would be more 
complicated than we have in Eq. (11/ 



III. EXAMPLES 

Our scope here is to apply the results obtained above 
to expand the Breit equation and the Salpeter equation 
to order 1/c 4 . For illustrative purposes, we consider the 
case of Coulomb particles with unequal masses and of 
charges t\ and £2, interacting through the potential 

1/(V . eie 2 eie 2 / (oi ■ r)(a 2 ■ r) \ 

V(r) = — — Jai • a 2 + - j , (13) 

where r = ri — r 2 and r = |r|. One should remember, 
however, that in general the complete form of the original 
interaction of two Coulomb particles also has to include 
many other components of lower order of magnitude in 
addition to the Breit interaction (I13[) . e.g., such as the 
intrinsic magnetic moment terms [8[ , the terms involving 
the electron self-energy and vacuum polarization (25j |. etc. 
Nevertheless, here for simplicity we restrict ourselves to 
the consideration of an expansion of the equations only 
with the potential in the form fp~3|) . 



A. The Breit Equation 

Let us consider the Breit equation Hip = Eip, which 
Hamiltonian is similar to the one in Eq. ([1]) if we put 



Using the expression (|6a|) , one performs the nonrelativis- 
tic expansion of the Breit equation to order 1/c 2 and, 
putting j3i = fa = 1, gets the Breit correction derived 
also in QED (see, e.g., [26|, [27j]). Based on Eq. (g]), which 
goes into (6a) in the case under consideration, an expan- 
sion of the Breit equation to order 1 /c 2 and the study of 
properties of its transformed Hamiltonian were carried 
out in [y, |8[ . Note that the Breit correction is divergent 
as it involves the Dirac <5-functions appearing because of 
the Coulomb potential in the original Hamiltonian and 
because of the same potential, the (^-functions, already 
together with their derivatives, will also appear and in 
the expansion terms in higher orders. 

Since the prescription © is just suitable for the case of 
the Breit equation, it is sufficient to work out the terms 
(|6bl ..i) with regard to the notations (fl4l . and thereby to 
get the l/c 4 -order part of the transformed Hamiltonian 
in an explicit form for this equation. In order not to miss 
any ((-functions and their derivatives in the final results, 
first one has to calculate these terms in momentum space 
and then to pass into coordinate space, as they usually 
do for the derivation of the Breit correction. Our main 
goal in this subsection is to obtain all the terms of order 
1/c 4 in the expansion of the Breit equation, which final 
forms we write below mainly in coordinate space. 

We start off with the terms ( |6jj ) 



16mfc 10 



(osy 



16TO2C 10 



(SO) 6 = 



16to^c 4 



(15) 



which provided that both particles are in positive energy 
states yield the correction of order 1/c 4 to the kinetic en- 
ergy. The rest of the members in ([6]) , except (|6a)) , lead to 
1 / c 4 -order corrections to the effective potential, form the 
interaction terms in the higher-order Breit-Fermi Hamil- 
tonian, and with respect to the contribution of (S S ) and 
(OO) can be divided into three groups. 

The terms from the first group give corrections to the 
Coulomb interaction. They are represented by (|6bl g, h) . 
For the first term in (j6b[) , after substituting the operators 
from (|14al b) and working out, one gets (h = 1) 



[(OS), (SS)}' 



(ei£2 



) 2 c 2 



(16) 



The calculation for the sum of the two terms in the first 
bracket of ((Sgl) gives the following formula: 



[(OS),[(OS),[(OS),[(OS),(SS)] 



(SS) 



£l£2 



(OS) = CQl -pi, (S O) = C£* 2 • P2, 



(14a) 
(14b) 



(OO) = — — I 0:1-0:2 H ^ )■ ( 14c ) 



-32[(OS) 3 , [(OS), (SS)} 



= 72e l£2 c 4 ir[pl,6(r)} + ~ 2ir(V6(r) x Pl ) • <7\ 



(rxpQ-^ .^(rxpi)-^ 



-Pi 



The even-even and odd-odd terms in this case denote the 
original interaction in H and commute with each other. 



15c 4 





D 2 — " 




p?, 


Pl! 






r 





(17) 
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We can work out the commutator as 

= -47rA<S(r) - 167iW<5(r) • pi + 4 





' 2 1" 




PI, 






r 





-Itt 



1 



— S(r)6V + _ U»J - 3^- 



P\p{- 



(18) 



The next two members in the second bracket in Q6g[ ) are symmetrical to the ones derived above and can be obtained 
from them by interchanging the indices 1 and 2. 
Then we have from the term (I6hl) 



[(0£),[(0£),[(£0),[(£0),(££)}]}] = / e 



47rei£2C 



(q 2 - Ha x ■ (q x pi)) (q 2 + 2ier 2 ■ (q x p 2 )) 



d 3 q 
(2tt) 3 



eie 2 c 4 ^ -47rA<5(r) - 87r(V<5(r) x p ± ) ■ a l + 87r(V<S(r) x p 2 ) • cr 2 



4tt 



1 



_ s(t)8» + s S* 3 - 3— 



(<ti x pi) 2 (er 2 x p 2 ) j 



(19) 



One can easily see that singular operators containing the 
Dirac (5-function and its derivatives enter into Eqs. (|17p 
and (HHJ). As already remarked, that is because the (££) 
term, in the original Hamiltonian, is represented by the 
Coulomb potential in the case considered here. Thus the 
terms (60 g, h) with allowance for Eqs. ([IB]). (|T7|) . $TS§ 
(and for their symmetric equations) give relativistic cor- 
rections to the Coulomb interaction in the transformed 
Hamiltonian. With neglect of the odd-odd operator (I14c[) 
in the Breit equation, only these terms form the higher- 



order transformed Hamiltonian, which can be applied to 
the equal-mass case as well. It should be noted that the 
term similar to (|6h| with (|19p was obtained for the effec- 
tive a e mc 2 Hamiltonian and discussed in [l9| . 

The second group of the terms, which is a little larger 
than the preceding one, comes from (|6c[ d, i) and thus 
is responsible for corrections that are conditioned by the 
Breit operator, namely, by (OO). 

The calculation of the second term from (I5c|) gives us 
such a result: 



[(0£),(00)]\ 



(eie 2 ) 2 c 2 



2 2 3 18 

13 - 3 ( r X Pi) ' CTl - 3( r X Pi) ' ^ - "3<Ti ■ T2 



7 4tt 



S(r) 



1 



2i 



(20) 



The terms of the next type from (|6dp . which are represented by the first and third ones, have a more awkward form 
and the calculations of them are rather laborious and lead to lengthy expressions. However, we can avoid their direct 
calculation and simplify the procedure instead, using the following relation (for the first term): 



[(00),[(0£),[(0£),(00)} 



= 2[(0£), (00)} 2 + + [(OS), [(OS), (OO) 2 



Here the first member standing on the right-hand side is already obtained in (|20[) . and at the same time the second 
member is much simpler than the one on the left-hand side. The calculation of it yields 



[(0£),[(0£),(OOf]\ = 



(eie 2 ) 2 c 2 



3 32tt 



6 



7 



-j + —z-8(r) + ^-(r x p x ) • (Ti + — (r x p x ) ■ er 2 + -j(Ti ■ cr 2 



1 4tt 



+ -7T*(r) 38 



2r l p{ + Zr 3 p\ 



'l u 2 



75V 



r Pi 



2S ij 



°"l CT 2Pl ~ -Pl)^! 



(21) 
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For the second type of the anticommutators from (|6d[) we have 



[(00),[(00),(OS) 2 } + } + 



(eie 2 ) 2 c 2 



x 3<5 lJ 



r 



1 647r r, x 1 / \ / ,13 

- + — <J(r) + -3 (r x Pl ) • (<n + <r 2 ) + — ^ ■ er 2 - 4 



1 4tt 



- j ff X 2 



, 12i 
Pi + — 



r Pi 



2i 



r • pj 



Art 



(22) 



The two members in the first bracket in ([61]), which are 
linear in (CO), e 2 terms, can be worked as follows: 

[(08), [(OS), [(OS), [(SO), (00)} + } + }+] + 



+8[(OS)\[(SO),(00) 

-eie 2 c 4 (2pi - (<x 2 ■< V ' 



r 



2p[ + («ri x V) 



+ (<7! X (2i P i - V)) 1 



(23) 



The rest of the members from (I6c[ d, i) are symmetrical 
to the calculated ones. 

Since the terms (|6bl g, h) involve the operators (I14a[ b) 
only, their calculation gives relatively simple equations, 
however, it is well seen that the terms calculated for the 
second group have lengthier expressions than the terms 
with Eqs. (flTj)) . fllTl) . and (H"9l) . The expressions obtained 
in Eqs. (|2T)|) - ([2"B"f have very cumbersome forms because 
the original terms (I6c[ d, i) contain six cr-matrices in our 
case and this leads to tedious calculations. 

At last, the third group of the expansion terms comes 
from (|6e[> . All of them involve both the even-even and 
odd-odd operators and therefore give contributions that 
are conditioned by both the static Coulomb interaction 
and the Breit operator. They are nonlinear in the original 
interaction, thus, are e 4 terms. 

The calculation of the first member in (|5e|) yields 

[(00),[(OS),[(SO),(SS)]]]+ 



2c 2 



d\d*k ^ 4 " £1£2 



(2tt) 6 J J H q 2 k 2 V k 2 

x {gV - (o-i x q)V 2 x qyy^+V' 

2(6i£ 2 ) 2 c 2 f 1 tor 3 

^ - |*(r)) (s^ - l£Wa} ■ (24) 



The expression for the second member has the form 
[[(SO),(SS)],[(OS),(00)]+] 



2(e ie2 ) 



1 87T 1 

-3 - — S(v) - -3 (r x Pl ) 



<T2 



5 lJ 



'l u 2 



(25) 



Notice that all the terms from this group are only spin- 
depended ones. So far as the e 6 term (JU is concerned, 
it equals to zero because (SS) and (OO) commute in the 
case under consideration. 

Thus, the terms [f6bl..j) with regard to Eqs. (JTSJ) - (|2"5|) 
form the l/c 4 -order part of the transformed Hamiltonian 
for the Breit equation. However, as shown in the previous 
section, one can change it a little bit, applying a unitary 
transformation with the generator in the form (Tj"2")) . Still, 
this transformation is not the only one that modifies the 
higher-order Hamiltonian. 

Consider, for example, a Hermitian even-even operator 
in the following form [28j ]: 



S, 



16c 6 



MOS) 2 , /3 2 (SO) 2 



.(££) 



(26) 



We bear in mind here that the operators (SS), (OS), and 
(SO) are defined by (|14a[ b). Under this assumption the 
function (l26l) reads 



Sep — 



16c 4 



fflP? , /?2Pl £l£2 



(27) 



Applying the procedure ([2} with the generating function 
in the form (|27[) to the transformed Hamiltonian for the 
Breit equation, and with neglect of the terms of sixth and 
higher orders, one obtains 



(iT-tr)' ~ #tr — ( —5 H % 

. raf 777,2 



8c 4 r 4 



32m 4 c 4 



Pl) 



Pl, 



£l£2 




1 




D 2 ^ 




r 




32m 4 c 4 


pI, 


21 r 





/3i/3 2 (t77 2 + 777 



32777^ 77720 



3„4 





D 2 






P2; 






r 





(28) 
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One can see that the term coming after H tr is similar, up 
to a sign, to the term obtained from (|6b| (see Eq. (|T6"])). 
and therefore this additional transformation removes the 
last one from H tr in the treated case. At the same time, 
given that the function ([26|i is taken with another factor 
such as — ib I (64c 6 ) instead of — i/(16c 6 ), one destroys the 
last member in (fT7|) . but the terms (I6b|) will be saved out, 
however, with a different numerical multiplier. 
Let us also consider two more generators: 



See — 
See = 



iftl 



16m 2 m 2 c 6 



[(0£),(00)} + ,(£0)}, (29a) 



[(£0),(00)] + ,(0£)}. (29b) 



Here the second operator is symmetrical to the first one. 
They can also be applied to the transformed Hamiltonian 
to modify its higher-order part. Using the function in the 
form of (I29bp . one obtains such new terms of order 1/c 4 
in the Hamiltonian: 



i[See,(££)} = 



Se 



Pi{C£f 



2m ic 2 



/3 2 (eie 2 ) 2 
8m 2 m 2 c 4 r 4 

-2(r x p x ) • a i + (r x p 2 ) • cr 1 



(30a) 



/3i/?2eie2 / T-,\\i 
3 j (f i x (2i Pl - V 



32m 3 m 2 c 

x (V 2 - (<t 2 x vy 



(30b) 



examples considered these procedures produce, as a rule, 
new terms instead of the removed ones. Note that these 
transformations only change the higher-order part of the 
approximate Hamiltonian. 

Obviously, using the same procedure with "even" gen- 
erators, one can easily modify the transformed one-body 
Dirac Hamiltonian ([7]) in a similar way. 



B. The Salpeter Equation 

Consider the Hamiltonian of the Hermitian part of the 
three-dimensional Bethe-Salpeter equation in coordinate 
space [12| ]: 



H = H 1 + H 2 + ± [(A+A+ - ArA 2 -) , V(v)]_ 



(31) 



where 



9 _i_ E a i H n 
H a = ca a -p a + /3 a m a c z , A* = — - 



2E a 



E a = (m 2 a c 4 + p 2 c 2 



1/2 



with a = 1,2, and V(r) is determined by (p~3]l . In order 
to reduce this Hamiltonian to the form ([T]), we evaluate 
the anticommutator, keeping the terms which contribute 
to order 1/c 4 in the expansion. Hence we have 



{££) = (ft+fr)^ 



+ 



8m 2 c 2 8m 2 c 2 ' r 
3/9iPi , 3/3 2 p| £ie 2 " 



32mfc 4 32m 4 c 



4„4 ' 



Se 



mo? 



2to 2 c 2 



32m 2 m 2 c 4 

x f 2 P J 2 - (cr 2 x vy 



(o-i x (2ipi - V)) 1 



p!, 



(30c) 



(Of) 
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4m 2 c ' 2r 

(aa ■ p 2 )p 2 ei£ 2 
8m 2 c 3 ' 2r 



*X 2 



Pi 


(ai • Pi)p 2 eie 2 " 


lC 


8m 3 c 3 ' 


r 










+" — ) 


+ 



We have in (I30b[) the expression that is opposite to the 
sign of the second term in the bracket in (l23|) . The trans- 
formation with the generating function (|29a[) leads to a 
similar modification of the transformed Hamiltonian but 
in its symmetrical part. 

Thus the higher-order Breit-Fermi Hamiltonian can be 
modified by the unitary transformations. We considered 
here these transformations with even-even functions that 
contain the operators determined by (IT4")) . but in general 
they may be other functions of r, pi, p 2 , cri, <t 2 as well. 
With due regard for transformations of this kind, one can 
destroy (or create) a number of higher-order terms in the 
transformed Hamiltonian. Meanwhile, as shown in the 



(£0) 



ca 2 • p 2 



4m 2 c 
ai-pi eie 2 



a 2 ■ P2 ("2 • P2)p! eie 2 
r 



8m^c 3 



4mic 2r 

(ai ■ Pi)p 2 ei£ 2 
8mfc 3 ' 2r 



l l«2 



(CO) = 



ftpj + fepj £l£2 
2r 



(5'-' 
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It is seen that to perform an expansion of the Salpeter 
equation, one has to use a prescription in a general form 
like (HJ, as the operators (OS) and (SO) do not commute 
in the case under consideration. Still, if we supplement 
expression ([6]) with the terms (lihl i), we will get a pre- 
scription for the transformation of the Hamiltonian (|3"Tj) 
to fourth order in 1/c. It should be noted, however, that 



as (OO) is of order 1/c 2 , the terms containing this oper- 
ator, except the ninth term in (|6aj) . make no contribution 
to the desired accuracy, and hence, can be ignored (all of 
the terms with a mass difference in the denominators are 
among them). Thus, after the evaluation, one obtains 
the transformed Hamiltonian from the Salpeter equation 
to order 1 /c 4 in the following form: 



„ _ R 2 , 0lV>i , fa , a x e l £ 2 01 + 02 
H tr w pmnc + + (/3i + 2 ) 



2mi 



Ar 16to 2 c 2 



(o-i • Pi), 



(fi • Pi), 



ei£2 



01+02 
16mim 2 c 2 

01 



(o-i ' Pi), 



32to?c 4 
Pi 

327711m 2 , c 4 
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' 768to 4 c 4 
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01+02 

192TOimfc 4 



i \ 6162 

(<Tl • pi), 

r 



ft 
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A Pi 

, 8mjc 2 



(32a) 



2?' 

(fi • pi), 
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(32b) 
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2r 
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(32c) 
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2r 









£ie 2 ~ 
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+ r J 



ftp? , ftpi 



2mi 2m 2 
(0-1 • Pi), 

symm. terms. 



ftpi + ftp| 

' 2m 1 2m 2 



(32d) 
(32e) 



Although this formula is only an intermediate result, the 
Hamiltonian written in such a form is quite transparent 
and suitable for a qualitative analysis. Besides, it is rep- 
resented in the form similar to the one in (j6|) with regard 
to the notations (TT4")) . therefore, is convenient to compare 
with the transformed Hamiltonian of the corresponding 
Breit equation. 

Indeed, expressing Eq. (|3"2")l in terms of the (SS), (OS), 
(SO), and (OO) operators from (fl4"l) . we can note that it 
is almost similar to equation ^ [29j]. The terms in (|3"2"1) . 



except the terms (|32d[) . are similar to the ones entering 
into the transformed Hamiltonian of the Breit equation. 
The e 4 terms, which we put together in (|32b[) and which 
are of order 1/m 3 , do not disappear when one particle 
is in a positive energy state and the other is in a nega- 
tive energy state. The similar terms together with their 
symmetrical ones also occur among the nonlinear terms 
(|6bl c, e). Note, all of them agree with those e 4 terms in 
the expansion of the Breit equation, which do not have 
a mass difference in the denominators, however, they are 
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four times higher. Meanwhile, the e 2 terms (|32c[) . which 
are of order 1/m 4 , exhibit another behavior. They have 
the factor /3x + /3 2 , and thus, become zero if one particle 
is in a positive energy state and the other is in a negative 
energy state. When both particles are in positive energy 
states, they coincide with the terms ( |6g] h, i) for the case 
of the Breit equation, but all of them are opposite in sign 
to those provided that the two particles are in negative 
energy states. As concerns the terms (|32d[) . there are no 



similar ones in Eq. (j6j) . One should remark that this part 
includes e 4 terms having the factor fix + (3 2 (the e 4 terms 
of (|32bj) are without this factor) and e 2 terms which do 
not disappear when the particles are in different energy 
states, as opposed to the terms (I32c|) . 

It is easy to check, however, that such terms can be re- 
moved by an additional unitary transformation like those 
considered in the previous subsection. Namely, applying 
the procedure ([2]) with the generator 



S = 



Pi 



p 2 . 



32c 4 
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r i r j \ 












r 2 j 




+- 



(33) 



to Eq. (|3"2"|) . one eliminates all the terms of (|32d[) together 
with their symmetrical terms. This transformation only 
destroys the terms (|32dj) . without changing the rest part 
of the approximate Hamiltonian and without producing 
any new terms to the desired accuracy. In fact, this gen- 
erator contains three functions similar to (|27l) and (|29|) . 
but with other factors. Here the first function cancels the 
first and third terms from (|32d[) with their symmetrical 



ones, and the next two functions cancel the rest. So, with 
neglect of the terms (I32d[) . H tI , in the case under consid- 
eration, involves only the terms similar in their structure 
to those occurring in the transformed Hamiltonian for 
the Breit equation. 

Using the results obtained above, we can thus represent 
the final form of the transformed Hamiltonian for the 
Salpeter equation, accurate to order 1/c 4 , as 



ftp? 



ei£2 



£i£2 
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-<7l • <T 2 
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Here the terms to second order are put together in (34a) , and the others form the 1 / c 4 -order part of the transformed 
Hamiltonian. We do not include the terms coming from (|32d[) in this equation and represent them apart: 



(Pi + AO 



(eie 2 ) 2 
32m?c 4 r 4 
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j-j (2pi + (<n x V) 1 ) (<r 2 x (2zp 2 + V)) J 



symm. terms. 



(35) 



The extension of the FW method applied to the Salpeter 
equation here gives the formula with these terms, which 
can be referred to extra terms. Yet, we cannot claim that 
they also arise in the nonrelativistic expansion obtained 
by a different method, e.g., by the method of expressing 
the small components of the spinor of the wave equation 
through its large components [1, Q . 

IV. SUMMARY 

In the present paper we have dealt with the extension 
of the FW method to a system of two spin 1/2 particles 
and, by means of which, we have briefly considered the 
problem of expansion of equal-time relativistic two-body 
wave equations to higher orders in 1/c. For the case of 
unequal masses, we have found the fourth-order part of 
the transformed Hamiltonian in a general form in terms 
of even-even, even-odd, odd-even, and odd-odd opera- 
tors. It turned out that in contrast with the case of the 
transformation to second order, the form of the higher- 
order transformed Hamiltonian depends on the order of 
application of the generators, i. e., as these functions do 
not commute with each other [30], certain extra terms 
can arise in the expansion of the original equation. Still, 
terms of this kind can be removed by an additional uni- 
tary transformation, which, however, does not change the 
transformed Hamiltonian to order 1/c 2 . The occurrence 
of extra terms in Ht r may be a feature of the expansion 
of relativistic equations to higher orders. Evidently, any 
of the model effective Hamiltonians giving a contribution 
to a 6 mc 2 to the energy levels of an arbitrary light atom, 
not only the approximate Hamiltonians that one obtains 
from the nonrelativistic reduction of the wave equations, 
is defined up to a unitary transformation with generating 
functions of order 1/c 4 . 

For illustration, we have considered the nonrelativistic 
expansion of the Breit equation and the Hermitian part of 
the three-dimensional Bethe-Salpeter equation with the 
Breit interaction up to the terms of fourth order. Due to 
the Casimir projection operators, the l/c 4 -order part of 
the transformed Hamiltonian obtained from the Salpeter 
equation is much simpler than the analogous part from 
the reduction of the corresponding Breit equation, and 



does not involve any terms with mass differences in the 
denominators or terms similar to those in (|6d[) by their 
structure. However, none of the e 4 terms from the ex- 
pansion of the Salpeter equation (|34)) disappears when 
one particle is in a positive energy state and the other 
is in a negative energy state, meanwhile all the e 2 terms 
become zero [3l| . Notice that a number of the 1 /c 4 -order 
expansion terms worked out for these equations agree up 
to a numerical factor, and in many features, are similar 
to the relativistic corrections in the a 6 mc 2 Hamiltonian 
for an arbitrary light atom (see, e.g., [l9|). 

In fact, the approach applied here to expand the wave 
equations is a straightforward way of deriving many 1 /c 4 - 
order terms giving oftmc 2 corrections to the energy levels 
for a hydrogen-like atom. We emphasize, however, that 
this work has been devoted mainly to an extension of the 
FW canonical transformation to the two-body problem. 
Although we have illustrated the application of the re- 
sults on examples of the Breit and Salpeter equations, 
the application to a real atomic system and calculation 
of the energy states are beyond the purpose of the paper. 
Because of the simple form of the original interaction we 
used in the equations, the obtained approximate Hamil- 
tonians are not complete and do not take into account 
many QED effects. The derivation of the total effective 
Hamiltonian contributing to a 6 mc 2 to the energy for a 
hydrogen-like atom by means of the considered method 
also requires the use of quantum field approaches and the 
diagram technique, which is of interest in itself and could 
be the subject of future investigations. 

The results of the work are applicable to any two-body 
system, besides, they can be generalized to the systems of 
three and more number of particles. They may be useful 
in the study of few-body systems when the application of 
quantum field theories causes difficulties and can also be 
of some interest from the theoretical and mathematical 
points of view in general. 
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